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Abstract 

The intermediate-state interaction and structure of amplitudes of complicated pro- 
cesses in medium (decays, reactions and the nn transitions) are studied. It is proposed 
to use the branching ratio of channels of free-space hadron-nucleon interaction as a test 
in the construction and verification of the models. The corresponding formulas for the 
processes in the medium are obtained. The connection between particle self-energy and 
amplitudes of subprocesses are analyzed as well. 
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In this paper the structure of amplitudes of complicated processes in nuclei is studied by 
the example of the decay in the medium. We briefly summarize the well-known results relating 
to the unitarity and analyze the value and role of intermediate particle self-energy. This point 
is important since an additional self-energy, as a rule, leads to suppress the process. The 
expressions for the branching ratio of decay channels in the low-density approximation are 
derived as well. The formulas obtained have a clear physical meaning (we emphasize this fact), 
which enables one to verify and correct the model. It is shown that the models of realistic 
processes in the nuclei should reproduce the branching ratio of channels of the corresponding 
free-space hadron-nucleon interactions. This can be considered as necessary condition for the 
correct model construction. 

Let us consider a free-space decay a — > n°n, for example, A — > ir°n. For a decay in the 
medium the annihilation 

a^n° + n^n° + M (1) 

and scattering 

a -> 7T° + n -> 7T° + n (2) 

channels take place. Here M are the annihilation mesons, n — > M and n — > n imply the 
annihilation and scattering of n in the medium, respectively. In the following the antineutron 
and a-particle are assumed non-relativistic and spinless. 

Let r o , T s and T t be the widths of the channels (1), (2) and the total width of the decay 
a — > 7i°n in the medium, respectively. For a decay in nuclear matter T s because the n 
annihilates in a time r a ~ 1/r, where V is the annihilation width of n in the medium. The T s 
is taken into account since the p-dependence of the results is considered. 

The interaction Hamiltonian is 

Hi = Hd + H, (3) 

H d = g*$Q* a + H.c, (4) 

where Hd is the Hamiltonian of the decay a — > n°n, H is the Hermitian Hamiltonian of the 
n-medium interaction taken in the general form. We focus on the intermediate-state interaction 
of the n and so the 7r°-medium interaction is inessential for us. 

The total decay width T t can be obtained from the unitarity condition: 

r* = ^(1- | Sa | 2 ) « ^2ImT u , (5) 

S = 1 +iT. Here To is the normalization time, To — > oo. The on-diagonal matrix element Tu is 
shown in Fig. la. It involves the full in medium antineutron propagator G m which should be 
calculated through the Hermitian Hamiltonian 7i: G m = G m (H). Since G m {7i) contains the 
annihilation loops, the realization of this scheme is extremely complicated. 
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Figure 1: a The on-diagonal matrix element Ta corresponding to the decay a — > 7r°n in the 
medium, b Same as a for the nn transition. 

In the optical model 

H -> (ReU n + HmUnM^n: (6) 

where U n is the optical potential of the n. Then G m (H) — > G m (U n )- If U n = const, 

G m (Un) = -- 1 — n ~ (7) 

(m is the antineutron mass) which greatly simplifies the calculation. 

However, in [1,2] it was shown that model (6) cannot be used for the calculation of r t by 
means of (5). This is because Eq. (5) follows from the unitarity condition and can be used 
only for the unitary ^-matrix, whereas Hamiltonian (6) is non-Hermitian; the corresponding S- 
matrix S = l + iT(Un) is essentially non-unitary. As a result, in the above-mentioned model the 
effect of antineutron absorption acts in the opposite (wrong) direction [1,2]. The importance 
of unitarity condition is well known [3-5]. Nevertheless, the model based on (5) and (6) is 
frequently used (see below) because it greatly simplifies the calculation. 

In this connection we continue consideration of some aspects of unitary models of multistep 
processes in the nuclear matter. The decay channels (1) and (2) were selected on the following 
reason. Let us consider the nn transition [6-8] in the medium followed by annihilation 

n^n^ M. (8) 

This is a simplest process involving intermediate-state interaction since the nn transition vertex 
corresponds to 2-tail diagram (see Fig. lb). As a result, Eqs. (5)-(7) lead to the simple formula 
r t = —2ImeG m e (e is off-diagonal mass [7,8]), as opposed to the decay a — > n°n shown in Fig. 
la. Besides, this process is of independent interest [9]. (We note that all standard calculations 
of process (8) are based on Eqs. (5) and (6) (see [8,9] for future references), which is wrong 
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on the reason given above.) Decays (1) and (2) are considered because they involve the same 
n-medium interaction as process (8). 

For the decay a — > 7r°n in the medium we calculate the T a , T t and the branching ratio of 
channels using the low-density approximation [10] for the n-medium interaction. The results 
obtained are generalized to the decay 

a^b + c (9) 
in the medium and reaction on a nucleon of the medium 

ai + N^b + c (10) 

followed by elastic and inelastic c-medium interactions. 

We calculate the width of the decay (1). Our plan is as follows. 1) First of all we construct 
and study the simplest process model suitable to the concrete calculations. 2) In the framework 
of this model we calculate the values listed above. 

In order to construct the model correctly, we consider first the free-space nN annihilation 
(see Fig. 2a) and the process on a free nucleon 

a + N^7v° + n + N^7i° + M, (11) 

shown in Fig. 2b (the free-space subprocess). The amplitude of free-space nN annihilation M a 
is defined as 

<M\ Texp(-i j dxHnN(x)) -l\nN>= N a (2n) 4 d 4 (p f - Pi)M a . (12) 

Here TinN is the Hamiltonian of the nN interaction, N a includes the normalization factors of 
the wave functions. M a involves all the nN interactions followed by annihilation including the 
nN rescattering in the initial state. 

We write the formulas corresponding to Fig. 2b. The interaction Hamiltonian is 

Hi = Hd + HnN- (13) 
In the lowest order in Hd the process amplitude is given by 

M 2b = gGM a , (14) 

^ = 7 \ 7 \o~i ' (1^) 

(Po — qo — m) — {p — qj 2 /2m + iO 

where p and q are the 4-momenta of the a-particle and 7r°, respectively, M a is given by (12). 

Since M a contains all the nN interactions followed by annihilation, the antineutron propagator 

G is bare. 
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Figure 2: a Free-space nN annihilation, b Free-space reaction a + N^ir + n + N^-K° + M. 

Consider now the decay (1) (see Fig. 3a). The background a-particle potential is included 
in the wave function of a-particle ^ a {x). The interaction Hamiltonian is given by (3). In the 
lowest order in TCd the process amplitude M 3a is uniquely determined by the Hamiltonian (3): 

M 3a = gGM™. (16) 

The amplitude of the n-medium annihilation M™ is given by 

<f\T exp(-i J dxH(x)) - 1 1 0n p ^ q >= N(27r) 4 5\p f - Pi )M™ (17) 

(compare with (12)). Here 1 0n p __ g > and < / | are the states of the medium containing the n with 
the 4-momentum p — q and annihilation products, respectively; iV includes the normalization 
factors of the wave functions. 




a b 

Figure 3: a Decay a — > it + n — > 7r° + M in nuclear matter. The antineutron annihilation is 
shown by a circle, b The same as a, but the antineutron propagator is dressed (see text). 
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The definition of the annihilation amplitude M™ through Eqs. (17) is natural. If the number 
of particles of medium is equal to unity, Eq. (17) goes into (12). The antineutron annihilation 
width T is expressed through M™ (see Eq. (24)). M™ involves all the n-medium interactions 
followed by annihilation including the antineutron rescattering in the initial state. Due to this, 
the propagator G is bare. As in (14), the antineutron self-energy S = since the interaction, 
which can generate S, is involved in M™. 

If the Hamiltonian H is expressed through the nN- and iVW-interactions, the amplitude 
M™ contains in medium niV-amplitudes and dressed propagators. In this case the following 
condition should be fulfilled: if p — > 0, the propagator is not dressed. However, our purpose 
is to study the general features of the simplest model suitable to the concrete calculations and 
not the medium effects. On this reason we consider the model (16) which contains the block 
M™ corresponding to the observable values. We would like to emphasize this paragraph. 

The fact that the propagator is bare is important and so we study this point in detail. Using 
the same Hamiltonian (3), we try to construct the model with the dressed propagator (see Fig. 
3b). Denote 

V = ReU n . (18) 

As in (6), the Un is the optical potential of n. We recall the Hamiltonian 7i involves all the 
n-medium interactions. In the Hamiltonian H we separate out the real potential V: 

H = V^n + n' (19) 

and include it in the antineutron Green function 

G d = G + GVG + ... = ^- v . (20) 

The amplitude M 3b of the process shown in Fig. 3b is 

M 3b = gG d M' a . (21) 

Since the amplitudes M 3a and M 3b correspond to one and the same Hamiltonian (3), M 3a = M 3b 
and 

G d M' a = GM™. (22) 

The propagator G d is dressed: E = V ^ 0. According to (22), the expressions for the propa- 
gator and vertex function are uniquely connected (if Hi is fixed). The "amplitude" M' a (V,H') 
should describe the annihilation. However, below is shown M' a and model (21) are unphysical. 
Comparing the left- and right-hand sides of (22), we see the following. 

(1) If the number of particles of medium n is equal to unity, model (21) does not describe 
the free-space process shown in Fig. 2b because Eq. (14) contains the bare propagator. 
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(2) The observable values (T, for example) are expressed through M™ and not M' a . Com- 
pared to M™, M' a is truncated because the portion of the Hamiltonian H is included in the Gj- 
M' a has not a physical meaning. 

(3) Amplitude (21) cannot be naturally obtained from the formal expansion of the T- 
operator Texp(-i / dx(V^ n ^n + H')). 

(The formal expression for the dressed propagator should contain the annihilation loops 
as well. In this case the statements given in pp. (1) and (2) are only enhanced. A particle 
self-energy should be considered in the context of the concrete problem. The dressed propagator 
arises naturally if V and H! are the principally different interactions and vertex function does 
not depend on V. In our problem one and the same field generates £ and M™.) 

(4) Equations (20) and (21) mean that the annihilation is turned on upon forming of the 
self-energy part S = V (after multiple rescattering of n). This is counter-intuitive since at the 
low energies [11,12] 

— > 0.7 (23) 

Of 

(at = a a + a s , a a and a s are the cross sections of free-space nN annihilation and scattering, 
respectively) and inverse picture takes place: in the first stage of n-medium interaction the 
annihilation occurs. 

The realistic competition between scattering and annihilation should be taken into account. 
Both scattering and annihilation vertices should occur on equal terms in M™ or G d . According 
to pp. (l)-(3) the latest possibility should be excluded. Model (16) is free from drawbacks 
given in pp. (l)-(3). It reproduces the ratio (23) as well. 
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Figure 4: On-diagonal matrix elements corresponding to the process (a — medium) — > tt° +(n — 
medium) — > (a — medium). 

Indeed, we calculate the width of the decay (1) shown in Fig. 3a. Put p = for simplicity. 
The pion wave function is <&(#) = (2g ^)~ 1 ^ 2 exp(— iqx). Using the amplitude (16), the decay 
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width T a is found to be 



r a = / d q —±— g 2 G 2 r, 

J 2q (2ny 



(24) 



where T is the annihilation rate of the n with the 4-momentum p — q. The corresponding 
amplitude M™ is given by (17). 

The number of exchanges with the medium or particle, which appears in (20), is very 
important but unobservable value. In particular, it is responsible for the particle self-energy 
and process suppression by the potential V. The low-density limit enables one to verify directly 
the condition (23). Indeed, in the low-density approximation T = vpa a and 



r a = J dqFi(7 a , 
g 2 G 2 vp 



Ti 



(25) 
(26) 



2g (27r)3' 

<?o = + m2 r> where m n is the pion mass. 

Let us calculate the total width T t of the decay a — > n°n in the medium. In the lowest order 
in Ttd the on-diagonal matrix element shown in Fig. 4 is given by 

dq i 



M, 



-I 



-GM?(0)Gg. 



(27) 



(27T) i q 2 -m 2 Tr + i0 

Here M™(0) is the forward scattering amplitude of n in the medium. We integrate over go and 
use the optical theorem in the left- and right-hand sides of (27): 



— 2ImM 3 = r t , 
-i-o 

l2ImM s m (0) = vpa t 
-i-o 

(T is the normalization time, T — > oo), resulting in 

T t = J dqF l( j t . 

Denoting |q|= q, one gets finally 



(28) 
(29) 

(30) 



r = 



In the differential form 



T a = JdgF(q)a a (q) 
JdqF(q)a t (q)' 

F=^-—. 

% 

dTg/dq = Oa 
dT t /dq o t 



(31) 
(32) 

(33) 
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where Eq. (23) was used. In a similar manner, one obtains 

= - < 0-3- (34) 
dT t /dq a t V ; 

For the description of the intermediate-state n-medium interaction in (25), (29) and (34) the 
low-density approximation has been used. 

We consider the free-space decay (9). Let for the decay in the medium the inelastic 

a^b+c^b+x (35) 

and elastic 

a^b + c^b + c (36) 

interactions of the c-particle take place. Here c — > x and c — > c imply the reaction induced 
by c-particle and scattering of c-particle in the medium, respectively. Let T r and T s be the 
widths of decays (35) and (36), respectively; ay and a s are the cross sections of free-space cN 
reaction and scattering, respectively; T t = T r + T s and a t = ay + a s ; q is the 4- momenta of the 
6-particle. With the replacement T a — > T r and a a — > a r , the relations (31)-(34) describe the 
decay channels (35) and (36). 

For the reaction (10) in the medium instead of (35) and (36) we consider the channels: 

a 1 + N^b + c^b + x (37) 

and 

(H + N^b + c^b + c. (38) 

In this case the similar relations take place as well. The corresponding changes in F are 
minimum and non-principal. For example, for the channels (37) and (38) we get 

da™/dq = 
dcr^/dq a s ' 

where a™ and a™ are the cross sections of the reactions (37) and (38) in the medium, respec- 
tively. 

Relations (33), (34) relate the in medium values dT a /dq, dT s /dq and dT t /dq with the free- 
space branching ratio of channels a a /a t and a s /a t . The similar relations (Eq. (39), for example) 
take place for the decays (35), (36) and reactions (37), (38) in the medium. The nuclear medium 
changes the amplitudes [13] and branching ratio of channels [14] in hadron-nucleon interactions. 
Besides, the strong antineutron absorption takes place. Due to this, at the nuclear densities 
T a /T t ps 1. However, if p — > 0, the relations (31)-(34) and (39) should be reproduced. 

The definition of annihilation amplitude through Eqs. (12) and (17) is natural since it 
corresponds to the observable values. In that event the intermediate particle propagator is 
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bare. In this sense the above-considered model is opposite to that shown in Fig. 1 which 
contains the full in medium propagator calculated through the Hermitian Hamiltonian. In the 
intermediate cases one can get a double counting or opposite error. As a rule, an additional 
self-energy leads to suppress the process and so we would like to call particular attention to 
this point as well as the unitarity (see paragraph below Eq. (7)). 

Finally, once the amplitudes are defined by (12) and (17), the propagator is bare; the 
processes on a free nucleon shown in Fig. 2 are reproduced; for the branching ratio of channels 
the relations (33), (34) and (39) take place. In fact, these relations should be fulfilled for any 
process model and can be considered as necessary condition for the correct model construction. 
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